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Abstract 

We present an integral representation to the quantum Knizhnik-Zamolodchikov equa- 
tion associated with twisted affine symmetry Uq{A2 ) for massless regime \q\ — 1. Upon 
specialization, it leads to a cojectural formula for the correlation function of the Izcrgin- 
Korepin model in massless regime \q\ — I. In a limiting case q —1, our conjectural 
formula reproduce the correlation function for the Izergin-Korepin model |^ at critical 
point q — —1. 

1 Introduction 

We shall consider the quantum Knizhnik-Zamolodchikov equation (gKZ equation) associated 
with twisted affine symmetry C/g(^2^'*) for massless regime \q\ = 1. In the earlier work |jl], ||, the 
qKZ equation associated with twisted affine symmetry Uq{A2 ) for massive regime — 1 < g < 0, 
was studied within the framework of the representation theory of Uq{A^^), Hou et al. 

gave free field realizations of the vertex operators, and realized integral representations of qKZ 
equation, as the trace of the vertex operators. Their integral representation gave the correlation 
function for the Izergin-Korepin model for massive regime —l<q<0. In this paper we 

(2) 

present an integral representation for qKZ equation associated with twisted algebra Uq{A2 ) 
for massless regime \q\ = 1. Let V = C^, and consider the i2-matrix R{(3) £ End(y V) 

(2) (2) 

associated with Uq{A2 ) (see sec. 2). The qKZ equation associated with C/<j(^2 ) \l\ = 1 is 
the following system of linear difference equations for an unknown function G2n{Pi, ■ ■ ■ ,132n) 
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taking value in the space V^"^^ . 

G2n{Pi, • • • , - Ai, Pj+l, • • • , f32N) 

= Rjj+iiPj - Pj+i - Xi)''^ ■ ■ ■ Rj,2N{f3j - I32N - Xi)~^ 

X Rij{(3i - pj) ■ ■ ■ Rj-i,j{(3j^i - (3j)G2N{f3i, ■ ■ ■ ,(3j-i,f3j,Pj+i, ■ ■ ■ ,(32n)- (1.1) 

Here Rj^kiP) £ End(y^^^) signifies the matrix acting as R{(3) on {j,k)-th. tensor components 
and as identity elsewhere, and diffence parameter A > 0. Upon specialization of difference 
parameter A = Stt and spectral parameters Pj = f3 + 3-Ki,Pj^]\f = /3, (1 < j < A^) , om' 
integral representation leads to a conjectural formula of the A^-point correlation functions of 
the Izergin-Korepin model in massless regime \q\ = 1. 

G2n{P + 3^i, •••,/? + 3^i, /3, • • • , (1.2) 

Actually, in a limiting case q ^ —1, our conjectural formula for the correlation function 
reproduce those of earlier work at critical point q = —1. 

In connection with "massless" ^KZ, we should mention about the work IH, @, |8|, in which 
the authors presented an integral representation to qWL equation associated with Uq{A^l^) 
for l^l = 1. In pioneering work [Q, F.Smirnov gave conjectural integral representations for 
the form factors of sine-Gordon model. M.Jimbo et al.'s integral representation Q gave a 
conjectural formula of the correlation function for the massless XXZ chain. The higher spin 
(spin 1) generalization of work [0, ^ was achieved in |^]. The [/^(^^^l^) generalization of work 
0, ^ was achieved in T.Miwa and Y.Takeyama |11] studied qWL equation associated with 



Uq{J^^\-^) for |g| = 1, and presented hyper geometric pairing in terminology of V.Tarasov and 



A.Varchenko llj. 



Now a few words about organization of this paper. In section 2 we formulate the system of 
difference equations. In section 3 we present an integral representation. In section 4 we show 
that our integral representation satisfies the system of diffence equations. In section 5 we give 
a supporting argument of our conjectural formula. 



2 Difference Equations 

In this section, we formulate the system of equations we are going to study, including the gKZ 
equation associated with f/g(^2 ) ~ 1- -'■^ ^^^^^ paper we parametrize the deformation 

parameter q of Uq{A2 ) by ^ as follows. 

q = -e-T, e>l. (2.1) 
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Consider the i?-matrix R{/3) G End(y V) of Ug{A2 ) acting on the tensor product of V = 

RiP)vj,^v,, = ^k,^Vk,R{P)i% (2-2) 

fcl,fc2 = ±l,0 

= -1^7i(«. (2.3) 

The scalar function k{P) will be specified below. It is chosen to ensure that the i?-matrix 
satisfies the unitarity (|2.6D and the crossing symmetry (2.7). Nonzero entries of the i?-matrix 
are given as follows. 

_ _ -e'^shf^) 

ii(/3)o,i:? = Rm% = ^ 

shf^(/3-27ri) 



— 1 — 10 -e2«sh(^^ 

sh(^(/3-27ri)j 
_ _ _ _ -shf^ 



sh (^{(3 - 2iTi) 



sh(i|)sh(^(/3-vri) 



sh(i(/3-27ri)) sh(i(/3-37ri) 



sh('^V-e"ich('^) +e^^ch('g 



sh ( - 27ri) ) sh ( ^(/3 - Svri) 



s h(f) (eich(f)-e^ch(g 
sh - 27ri)) sh - 3^i) 



sh - 27ri)) sh - Svri) 



-/3 + 27ri 



sh ( - 27ri) j sh f - Svri) 



sh(^) sh(^)sh(^ 

^{PJOfl TZ NT TZ x TZ T- (2.4) 

sh f ^ (/3 - 27ri) j sh f ^ (/3 - 27ri) j sh f i - Sttz) J 
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Our ii-matrix R{(3) gives the Boltzmann weight of the Izergin-Korepin model The i?-matrix 
R{(3) satisfies the Yang-Baxter equation. 



RuiPl - (32)Rl3{f3l - f33)R23{(32 - P3) = R23{fi2 " p3)Rl3[Pl " /33)i?12(/3l " P2): (2.5) 

the unitarity condition, 

Ri2{m2i{-fi) = id, (2.6) 

and the crossing symmetry, 

Rifi+'^t'"' = (ie-«)'-i"xR('-,S+^V"*. (2.7) 



27,, V 2 

Let A'' be a non-negative integer. Consider a F^^^-valued function G2n{Pi, ■ ■ ■ , I32n), depend- 
ing on the 'spectral parameter' • • • , (32N- 

G'2Af(/3l, • • • ,/32Ar) = ^ Vj^i^---i^Vj2r^G2NiPl,--- ,P2N)ji,-,j2N- (2-8) 

jl,--- J2JV=±1,0 

We study the following system of difference equations for G2n{/3i, ■ ■ ■ 1 P2N)ji,-- ,j2N involving 
parameter A. Hereafter we assume that 

2(^ + l)7r > A > 0, C>1. (2.9) 

i?-matrix Symmetry 

G2n{(3i,--- ,(^s+1,Ps,--- ,/32N)ji,-,js+l,js,---,j2N 

Yl R{Ps - Ps+iy/JjJ^_lG2N{Pl, - ■ ■ ,Ps,Ps+l,--- '/32Af)ji,...,j^J^_^^,... j2jv-(2-10) 

Cyclicity Condition 

G2n{Pi,--- ,P2N-l,P2N — i^)ji,-- ,32N = G2N {P2N , P2, ' ' ' ^ P2N -l)j2N Ji,--- ,j2N-l- (2-11) 

Recursion Relation 

G2n{Pi, ■ ■ ■ ,Ps-i,P,l3 + 37ri,/3s+i, • • • , P2N)ji,--- ,js-i,j-j,js+2,- J2N 

= Gj G2N-2iPl, - ■ ■ ,Ps-l,Ps+l,--- ,(32N)ji,-,js-i,js+2,--,j2N- (2-12) 

Here Gj are specified below. See ( [4.27| ), ( [4.28| ) and ( [4.29| ). Note that the above system of 
equations involve only the functions G2n{Pi, • ■ ■ ^ (32N)ji,--- ,j2N with fixed value of 'spin' ji + 
\-j2N = 0. 
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3 Integral Representation 

The purpose of this section is to present an integral representation. In what follows we 
use the Multiple-Gamma functions Tr{x\uJi ■ ■ ■ uJr) and the Multiple- Trigonometric functions 
Sr{x\uJi ■ ■ ■ UJr) introduced in [|7| as follows. 

(-l)*" f e^^Hog(-t) dt 

logTr{x\uJi---UJr) =J — Br,rix\uJi, ■ ■ ■ , UJr) + ^ y , (ReX > 0), (3.1) 

r! Jc llj=i(l - e '^J*) 2mt 

where 7=Euler's constant, the integral contour C is shown in below Figure (Contour C), and 
Multiple-Bernoulli polynomials Br^r{x\uJi ■ ■ ■ uJr) are given by 

T-rr / =y2—,Br,n{x\uJi,--- ,^^r)- (3.2) 

The Multiple-Gamma function Tr{x\uJi ■ ■ ■ ujr) is an entire function of x. Tr{x\uJi ■ ■ ■ uJr) is 
meromorphic function with poles at x = niwi + • • • + rirUJr, (n-i, • • • , n,. < 0). 
Let us set Multiple- Trigonometric function by 

Sr{x\uJi, • • • , UJr) = Tr{x\uJi, • • • , UJr)~^Tr{uJi + ■ ■ ■ + UJr — x\uJi, • • • , UJr)^~^^ ■ (3-3) 



Properties of these functions are listed in Q. For examples they enjoy 

Trix + UJr\uJi, ■ ■ ■ ,UJr) 1 

rj.(x|a;i, • • • ,UJr) Tr-l{x\uJi, ■ ■ ■ ,UJr-l)^ 
Sr{x + UJr\uJl, ■ ■ ■ ,UJr) 1 

Sr{x\uJi, ■ ■ ■ ,UJr) Sr-l{x\u}i, ■ ■ ■ jW^-l)' 

and 



(3.4) 
(3.5) 



^ / , N ^ lT(x/uj) „ , , , /"KX\ 

Ti(xuj)=uj^-^ A . ^i xw =2sin — . 3.6) 
V27r \ uj / 



Here T{x) is Gamma function. As x — > oo (itlmx > 0), the function S2{x\ujiuj2) behaves as 
follows. 

log5,(x|a.,.,) = ±v^^ f-^ - ^^x " f ^ + ^ + ) + o(l). (3.7) 

\2u!lUJ2 2uJlUJ2 12 \UJ2 UJl J J 




Contour C 
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3.1 Auxiliary functions 

The integral formula involves certain special functions k{(3), p{P),ip{P),ijj{(3). 

• k{P) : Let us set 

'g2(-^^ + 27r)S2{-i(3 + 37^)g2(^/? + 57r)S2{iP + 67r) , , c ^ ,0 . « ^^ 
"^^^ = S,{iP + 2.)5,(i/3 + 3n)S,{-iP + 57r)g,(-i/3 + 67r) ' ^^'^^^ = ^2(-|2-e, Ott)). (3.8) 

The function K{a) satisfies the following difference equations, which ensure the unitarity and 

the crossing symmetry of i?-matrix. 

k{P)k{-P) = 1, (3.9) 

sh(^(/3-7ri))sh(|) 

k{P)k{/3 - Siri) = — ^ 4 — ^ ^. (3.10) 

sh (^{(3 - 2m)] sh (^{P - SmU 

• p{f3) : Let us set 

... ^ S^{-iP + 2iT)S^{-ip + 37r)53(z^ + 37r + \)S^{if5 + 27r + A) 

^^^'^ S3(-z/? + 27re)-S3(-z/? + 7r + 27r053(i/3 + 7r + 2^e + A)53(z/? + 27re + A)' ^ ' ' 

{S^{x) = S^{x\2<AM))- 

The function p{a) satisfies the following difference equations. 



• (^(a) : Let us set 



= «(/?), p{iX-(5)=piP). (3.12) 



'^^"^ 52(ia + 7r|A,27r052(-ia + 7r|A,27rO' ^^'"^^^ 
The function ip{a) satisfies following difference equations. 

ifia) = ^i-a), (3.14) 



ip{a — iX) 
tf{a) 



sh 



(^(a-TTi)) 



sh ^^(a + TTi — Xi) 



ipia±27r^i) _ sh(^(a±7rz)) 



(3.15) 
(3.16) 



<^(") shQ(aT7^^±27re^))' 
The function ip{a) has poles at 

(^(a) = ibi (niA + n227r^ + tt) , ni,n2>0. (3-17) 

The function (/'(a) is evaluated as follows. 

|</?(a)| < Const. e^p( 7—^- ^|q:| ] , |a| — >■ 00, (3.18) 

V 2^X J 

a/^ 1 

V 27r J- 

(p(a) = ^ ,„ , , „ — -r X h • • • , a — TTZ. (3.19) 

^ iS'2(27r|A,27rO a - ttz ' ^ ^ 
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tp{a) : Let us set 



The function "0(0!) satisfies the fohowing difference equations. 



i^ia) = (3.21) 
V'(a-iA) _ sh(^(a + 27ri)) 



^(") sh(^(a-Az-27ri)y 

TT 

a' 



^P{a±2n^i) shQ(aT27ri) 



^(/^) sh(^(a±27ri±27rei))' 



and 



(3.22) 
(3.23) 



(/j(a + TTi)ip{a)ip{a — vri) (3.24) 
= 2-» {sh (i(a + 2.i)) Bh (ii) .h + 2.,:)) }"' 

X |sh(^^(a-27ri))sh(^^a)sh(^^(a + 7ri))}"\ (3.25) 
The function V'(ck) has poles at 

a = ±i(niA + n227r^-27r), ni,n2>0. (3.26) 
The function ip{a) is evaluated as follows. 

|V'(a)| < Const, exp f — ^^ "*"o'^\^ ^^'^ 1*^1 ) ' 1*^1 ~^ (3.27) 
3.2 Integral Representation 

In this section we present the integral representation for G2n{Pi, ' " ^ P2N)ji,- ,j2N^ which is 
our main result of this paper. At first we demonstrate integral representations for N = 1 case, 
G2(/3i,/?2)o,o,G2(/3i,/32)i,-i, and G2(/3i,/?2)-i,i- 

C.(A.A,,^2Vf»h(|)p(A-A)/^|l/^|| 

¥'(«!- a2)V'(ai - a2) J] V^(«s - A) > e^^"^+"'"^^"^'^/iA(ai - "2) (3.28) 

s,t=l,2 J 



X < 



X sh (^^(«i - "2 + 27ri)^ sh (^^(""i + "2 + ttz)^ sh (^^(«2 - A + 7ri)^ sh (^^(-«i + P2 + vri) 
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X < ^{ai - a2)V'(ai - "2) ^^^^ ~ l^t)\ ^ ^^^^hx{ai - a2) (3.29) 

s,t=l,2 J 

X sh - «2 + 27ri) j sh + "2 + vr?) j sh " A + ^0 j sh ^^(o;2 - /?! + vri^ 



and 



X < (/p(ai - a2)'i/'(«i - «2) (^(os - > e"2e/5i/j^(Q,^ _ 0,2) (3.30) 

s,t=l,2 J 

X sh (^^(«i - "2 + 27ri)^ sh ^^(-«i + "2 + vri)^ sh ^^(""i + P2 + Tri) j sh ^^(-"2 + /?2 + vrf 
X (^e^"'sh(^^(a2-/3i+7rz)^ - e^"^sh (^^(/3i - ai + ttz)^ ^ , 
Here we have set 

hx{a) = sh Q(a - 2m)j sh Qa) sh Q(a + 27ri)) . (3.31) 

The integral contour C is taken along a path going from —00 to 00 in such a way that — tt < 
Im(ajk — Pj) < TT for all k,j. 

Next we present the integral representation for G2iv(/?i,-" ^ P2N)ji,--- ,j2N- This is Main 
Result of this paper. Given a set of indices ji, • • • ,j2N G {— 1,0, 1}. Let us set index set 

Aj, (j = -1,0,1) by 

Aj = {l<s< 2N\js > j}. (3.32) 

The index set A-i = {!,••• ,2A'"}. Let us set a_i,s = /?s, (1 < •§ < 2A'"). Let us define the 
kernel function ^'({a}|{/3}) by 

27V 

*(WI{/?}) = n n n n vi^o,s-ai,tMao,s-a,,t) 

s=ij=o,iteAj seAo,teAi 

^ n n '^K> - aj,t)V'(«j,3 - "j- *)• (3.33) 



j=0,l s,teAj 
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Here we have use the function {p(a),'ip{a) in (3.13), (3.20). Let us define the auxihary functions 
9fc({a,} •=-!), (A; = ±1,0) by 



5r_i({a_i}) = 1, 



5o({a-i,ao}) 



sh ( ^(ao — + 



(3.34) 
(3.35) 



gi{{a-i,ao,ai}) 



ie"Tsh(^ 



sh ( ^(ao — oi-i + vri) ) sh ( ^(ai — a_i + vri 



X ^e2«°°sh (^^(ai - «-i + ^^)) - e^s^^sh (^^("-i - ao + ni) ] ] . (3.36) 
Let us set the integral representation G2N by 



G2n{Pi,--- ,/52iv)ji,-- J2JV 



2Ar is 



2muN-,) -Q p(;3.-A)nn 



l<s<t<2iV 



s=l i=0 



da 



27ri 



^({a}|{/?}) 



s<t s<t s>t 

n *.({aj,s}f=_i) n n (^^"^'^ ~ ^ ^^^^ (^^""^'^ ^ ^ 
n 

aeAQ,teAi ^ ^ 

n sh f (ai,4 - ao,s + 27ri) j sh ^ 

fj [] [] sh (jii-Ps + aj^t + TTi)] H sh + /?s + 

s=l 7=0,1 te.A, ^ ^ ^ tsA. ^ ^ 



ai,t + 27ri) ) sh ( — (-ao,s + ai,t + vri^ 



,2? 



-ai J + Qo,s + T^i, 



(3.37) 



Here we have used p{f3) in ( |3lID , ^'({a}|{/3}) in (Q), /iA(a) in (|OTJ), ({a^- J^4_;^) in 
( 3.34| ), ( p. 35 ), ( 3.361 ). Here the integral contour C is taken along a path going from — oo to oo 
in such a way that — vr < Im(afc — /3j) < vr for all A;, j. We understand Go = 1. 



Let us check the convergence of the integral of G2n- Let us set 
Jx{a) = v3(a)sh ^-^(-a + 7ri)^ . 



(3.38) 
(3.39) 
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Let us consider the evaluation of |afc,t| — * oo. In our integral representation, factors involving 
the variable a^^t consist of following three parts. 

5it({ai,t}f=_i) 

X llJx{-f3s + ak,t)llJx{-ak,t + l3s) II (3.40) 

s<t t<s j=0,l 

s<t i=0,l t<s j^k 

For \a\ oo, the factors Ix{a), Jx{a) are evaluated as follows. 

Ix{a) = Const. exp (^^(^ - l)|a|^ , (3.41) 
J^{a) = Const, exp ^^(1 - ^)|a|^ . (3.42) 

We have fi'jt ({oj,t}jL-i) ~^ ^' i'^k,t oo). Let us consider the evaluation for variable |afc,t| — *■ 
oo in the integrand of our integral formula. Because of spin condition, ji + • • • + j2N = 0, 
there exist 2N factors of type ( 3.3g| ), (2A^ — 1) factors of type ( |3.38| ), and one S'jt ({oj,t}jLi)) 



in the integrand. Using above estimates ( p.4lD ( 3.42| ), we know convergence of the integral 
representation upon condition ^ > 1. 



4 Proof of Difference Equations 



In this section we prove the integral formula ( 3.37 ) satisfies the system of difference equations 
(|2l0|) , (HID, and (ilD. 



4.1 i?-matrix Symmetry 



In this section we prove iZ-matrix symmetry (2.10). Let us set 



G'2Ar(/3l, • • • ,/32Af)ji,- J2JV = Y\. P((^s-Pt) ^G2NiPir-- ,P2N)ji,-,j2N- (^-l) 



l<s<t<2Af 



Because of the properties of = n{P),p{iX — (3) = p{P), the equation ( 2.10| ) is reduced to 

the same equation for G2N where R is replaced to R. 



Jl,--- Js+l Js,--- J2JV 



Yl Rif^s - Ps+iyj'Jj^llG2N{Pl, - ■ ■ ,Ps,Ps+l,--- >/327v)ji,...jVJ^^^,...j2^. (4.2) 
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• N = 1 Case. 

First we demonstrate how to prove ( |2.10| ) in simple cases, = 1. There exist 3 cases to 
consider. We prove ( 2.10| ) by checking every cases. Let us consider the case, 



G2(/32,A)-i,i. (4.3) 
Considering the integrand of (LHS)-(RHS), we get 

xhxiai -02) Int(aia2|/5i/32)-i,i, (4.4) 



where we set 



Int(aia2|/3i/52)~i,i (4.5) 
sh ( — (ai — Q2 + 27ri) ] sh ( — (— ai + 02 + vri) 



2^' 7 \2(, 

X jish (^(-ai + /3i + sh (^(~"2 + + 

X (^e^^^sh (^("2 -P2 + vri)^ - e^"^sh (^(/32 - ai + vri) 

- Ripi - /32)l;:}ish (^(-ai +f32 + TTi)^ sh (^(-"2 + /?2 + vri)^ e"^^^ 
X (^e^"^sh (^(02 vri)^ - e^"'sh (^(/^i - "i + vri) 

- i2(/3i - /32)?;° ish (^(a2 - /3i + ^i)) sh (^(-"i + P2 + vri)) ei^^^+^^-^i"^^) 

- RiPi - P2)i'Aisii (^(ai - /3i + TTi)^ sh (^^(a2 - A + vri)^ e"^^^ 

X (^e^-ish(^^(a2-/32 + vri)) -e^"'sh(^^(/32-ai+7ri)^)|. (4.6) 

Because of antisymmetric relation hx{a) = —hx(—a), we get the equation ( [4. 3D from the 
following relation of trigonometric function. 

Int(aia2|/3i/32)-i,i - Int(a2ai|/3i/?2)-i,i = 0. (4.7) 
Other cases are shown as the same manner as the above case. We compare integrands directly. 



N >2, General Case 
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Next we prove general case. There exists the fohowing decomposition of the integral represen- 
tation G2Ar(/3l • • • (32N)ji-j2N- 



J{a[ ■ ■ ■ a'j^^j^^^^2\Ps, Ps+i)js,js+i 



X Sym(- • • , a'l • • • , /3s, • • • )• 

Here J{a[ ■ • • a^-^+^-^^^+al/^s, is the integrand of G2{f3s, Ps+i)js,js+i^ 

Here G2(/3i , P2)js,js+i is formally introduced by using the equation ( f3.37 ) for not only js+js+i '■ 



but also js + js+i = ±1, ±2. The function Sym(- • • , • • • . . . | . . . , /3^, /3^+i, • • • ) 

is symmetric with respect to a'^ <-> a[ and /?s <-> /^s+i. Therefore i?-matrix symmetry ( [4 .21) is 
reduced to the following relations of trigonometric functions. 

Yl -''Ka)'--- '«ao.+j,+i+2)l^s+i,/3sk+i,i. (4.10) 

The above equations for js + js+i = 0, have appeared in = 1 case. For js + js+i = ±1, ±2, 
we have to do tedious checking. 

4.2 Cyclicity Condition 

Because of the relation p{l3) = p{i\ — /?), cyclicity condition ( 2.11| ) is equivalent to 



G2Ar(/3l,--- ,/32Ar-l,/32Ar - ,j2N = G2NiP2N , Pi, ■ ' ' y P2N-l)j2N,ji,- j2N-i- (4-11) 

N = 1 Case 



First we demonstarate how to prove (2.11) in simple case, N = 1. Let us consider the case, 



G2(/?i,/32 - = G2(/32,/3i)i,-i. (4.12) 

In the RHS we change the variable ai — > ai — iX, a2 ^ CX2 — iX. The integrands of the LHS 
and the RHS coincides because of the relation. 



(f{a-iX) ^'^\2( 



sh ( ^(a — m) 



sh(^{a + m- Xi) 



(4.13) 
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We have to check that the contours for the LHS and the RHS are the same. Consider the 
LHS. The points aj = /3i + vri which appear in poles of the kernel (/?(aj — are actualy not 
poles, because there exist zeros from the factor sh ^^(— ai + /3i + vri)^ sh ^^(—02 + /3i + tt*)) • 
Therefore the integral J^j^.^dai J^,^.^da2 can be deformed to J^,dai J^da2- The parameter 
condition, 

A < 2^(^ + 1), (4.14) 

ensures that there is no poles in deforming strips of the integral. Therefore there exists no 
difference in integral contours between the LHS and the RHS. Other cases are proved as the 
same manner as this caes. 



• General N > 2 Case 

Let us consider general case ( [4.11 ). For simplicity we consider the case j2N = +1- We make 



the following change of integration variables : 

ao,2Af ao,2N — ^A, ai^2N ai,2Ar — in the LHS, (4.15) 

and 

in the RHS. 

After changing variables, the difference between the LHS and the RHS appears only in the 
factors involving variables /32Af; cto,2Af; CKi,2Af- The factor of LHS, which involves the variables 

P2N, ao,2Af, ai,2Af is 

2N-1 js 

5'j2Jv({"j,2Ar}i=±l,o) Y\. n '^^("i''^ ~ ^2Ar + «A) J| JxiP2N - aj^2N) h{ao,2N - ai^2N) 

s=l j=0 j=0,+l 
2N-1 2N~1 js 

X n n MPs-aj^2N + iX) n n n h{aj,s-ay^2N + i>^), (4.17) 

s=l j=0,+l s=l jr-=Oj'=0,+l 

Those of RHS is 

2N-1 js 

5i2Jv({"j,2Af}j=±l,o) 

n llMP2N-aj,s) n '^^(^2iV - aj,2Af) h{ao,2N - ai,2Af) 
s=l j=0 j=0,+l 
2N-1 2N-1 js 

X n n Mo^j,2N-Ps) n n n ^A(aj',27v-aj>). (4.18) 

s=l j=0,+l s=l j=Oj'=0,+l 
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Here we have used Ix{a), J\{a) in ( 3.38| ), ( p.39| ). Using the following property, 



Ix{a) = Ix{i\-a), Jx{a) = JxiiX - a), (4.19) 



we know that ( [4.17| ) and ( [4.180 are the same. We have shown that the integrands of the 



LHS and the RHS are the same. We also have to check that the contours for the LHS and 
the RHS are the same. Consider the LHS. The points aj^2N = Ps + T^hiJ = 0, +l;s = 
1, • • • , 2N — 1) which appear in poles of the kernel <f{aj,2N — Ps) ai^e actualy not poles, because 
there exist zeros from the factor Yll=i^ Y\j=Q ^ish (^■^{—aj^2N + Ps + "^i)^ ■ Therefore the 
integral Jc+ix^'^o,2N f(j_^_ix^'^^:'^^ can be deformed to fQdao 2N fcd'0ii.2N- The parameter 
condition A < 27r(,^ + 1) ensures that there is no poles in deforming strips of the integral. We 
have proved the equation ( [4. Ill ) for j2N = +1- For j2N = 0,-1 case we can show the relation 



( 4.11 ) as the same manner as the case j2N = +1- 



4.3 Recursion Relation 

It can be shown that Recursion relation ( |2.12| ) are consequence of i?-matrix symmetry ( p. 10)) , 
unitarity crossing symmetry (|2.7| ), and 



G'2Ar (/?i , • • • , /32Af-2, /?,/? + 37ri)jj^...j2jv_2,j,-i = CjG2N-2iPi, - ■ ■ ^ P2N -2) n,- ,j2N -2' (4.20) 

We shall prove ([4. 201) . Consider the integral formula (|3.37|) . The factor p{P2N-i — P2n) has 
zero at P2N = P2N-1 + Svri. We see that 

-1 



I / /IT 



p(/?2Af-i - /?27v) = p(37ri)^ (^cos I^^J sin J (/32Ar - /?27V-i - Svri) H . (4.21) 

The integral may have a pole at P2N = P2N-1 + Svri because the integral contour (—00, 00) is 
pinched by the pole of the kernel function ^'({a}|{/3}). We explain this procedure using simple 
example. For regular function f{ai,a2), the following estimate holds. 

dai f da2 ^ 



C-^ 27ri Jc^ 27ri "2) - ai + 7ri)(ai - 02 + 7ri)(a2 - P2 + T^i) 

= —m ^ , o -N + Regular, (/?2 ^ /3i + Svri). (4.22) 

(Pi - P2 + 3vrz) 

Here the contour Ci is taken along a path from — cxd to 00 in such a way that — vr < Im(a2 — /?2), 
Im(ai — a2) < vr, and the contour C2 is taken along a path from —00 to cxd in such a way that 
— vr < Im(Qi — 02), Ini(/3i — ai) < it. We will check this procedure for our considering case, 
and calculate the residue of G2N ■ 
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• iV = 1 Case 

First we demonstrate simple case, N = 1. 

G2{Pi,/32)j,-j, /32^Pi + Siri. (4.23) 

The integral J^dai is pinched by the pole of f{ai — and the integral f^da2 is pinched 
by the pole of f{a2 — ai). Because S2{uJi\toiUJ2) = ^ The kernel function (/^(ai — a2)V'(cti ~ 
a2) Wg 2 '^{'^s — Pt) behaves for /?i — > ai + iri, ai — > a2 + 27ri, a2 ^ P2 + "^i- 



- 3 

K 



ip{2mfilj{m) J ^^ I Z5 ^7 ^ ^7 o ^ + " " " • (4-24) 



1 



zS'2(27r|A, 27r^) I - ai + 7rz)(Q;i - 0:2 + 7rz)(Q;2 - /52 + ttz) 



Note that 



^(27rz)V(^0^A(7rz) = (2V(0)sh2 sh ' . (4.25) 

We know the following is independent of the spectral parameter /?. 

G2{P,P + 3Tri)j,.j = Cj, (i = ±l,0), (4.26) 



where 



V2e~ 


772 






V2e" 




V2e~ 







-3 

_ v-^ ^ /--s .2 /^7ri\ p(37rz) ip{27Ti) 

- V UJ 52(27r|A,2.0^ ^(0) ' ^'-''^ 

-3 

_ V - 1'^: . fTri\ pi^^-i.) y(2^i) 

- ^^^/^^'UJ 52(2.|A,2.0^ ^(0) ' ^'-''^ 

.3 

_ V- " /'-s ,1^2/^^ A P(37r0 y(27ri) 



• General N >2 Case 

We show that general N > 2 case is reduced to iV = 1 case. Consider the integrand of the 
following integral formula. 

G2n{Pi,--- , P2N-2, P'l, l3'2)ji-j2N-2,j-j- (4-30) 
Let us set integral variables, 

{a'l,^} = {{ak,2N-i}i=odoik,2N}kio}- 
The factor interacting both aj^g, (s = 1, • • • , 2N — 2; j = ±1, 0) and I3[,P'2, a[, 0I2 is given by 

n p{i5s - 0x)p{iis - 02)JxWs - oc'iVxWs - a'2) 

s=l 
2N-2 js 

^ n llM(^j,s-Pi)M»j,s-P2)h{aj,s-a[)Ix{aj,s-a'2). (4.31) 

s=l j=0 
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The above factor ( [4.31 ) simplifies for /?2 ^ 02 + T^h (^2 ~* + '^1 ~^ + '^h 

(4.31)^2-28(^-1), (4.32) 

because of the following relations. 

p{(3)p{P - 37Ti)Jx{P - vri) Ja(/3 - 2^i) = 2-^, (4.33) 

Jx{P)Jx{P - 37ri)/A(/3 - 7ri)h{/3 - 27n) = 2^^^ (4.34) 

This simplification teaches us that general N > 2 case is reduced to the simplest = 1 case. 



We have proved recursion relation ( 2.12 ) 



5 Supporting Argument 

The purpose of this section is to give a supporting argument on physical meaning of our 

integral representation. Let us consider the integral formula for A = Sir. In this case following 

simplifications occure. 

^ S2{-if3 + 27r\6n, 2nC)S2{-iP + 3^|67r, 27rO 
^^^^ 52(-i/3 + 2<|67r,27r052(-i/3 + 2< + 7r|6^,2O' ^ ' ' 



tpia) = f{a) X ^ 1- ^ ^. (5.2) 

4sh 2? (-a - 2Tri) )sh{^{a- 2-111) 



and 



2gV - ^-.vy--\^25> 

For the case A = 37r, our integral formula ( |3.37| ) gives a conjectural formula for the Izergin- 
Korepin model Hou et al.|]l|, ^ considered the correlation functions for the Izergin-Korepin 

model for massive regime — 1 < g < 0, within the framework of representation theory of 

(2) 

Uq{A2 ). They gave free field realizations of the vertex operators, and realized integral repre- 
sentation of the correlation functions, as the trace of the vertex operators. In the limiting case, 
our integral representation reproduce the correlation function for the Izergin-Korepin model at 
critical point q = —1, which was derived by Hou et al.[|^, ^. This gives a supporting argument 
that our integral representation (3.37) gives a conjectural formula of the correlation function 



of the Izergin-Korepin model at massless regime \q\ = 1. 

In the special case where A = 37r and ^ co, the auxiliary functions gk{{c(j}j=-i), {k = 
±1, 0) tend to gk{{»3}j=-i), (k = ±1, 0) given by 

9-i{{a-i}) = 1, (5.3) 

5o({a-i,ao}) = : (5.4) 

ao — ct-i + TTi 

~ ,r 27r(2a_i - ao - ai) 

fifi(|a_i,ao,ai|j = 7 : ^7 : r- (5.5) 

(ao ~ CK-i + vrzj(ai — a_i -|- m) 
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The intergal representation G2N, ( 3.37] ) becomes as follows. 

n 



l<s<t<2N i ^ Gtt ^ ^) ^ \ Gtt ^&)^ \ &^ 2 j M ^ ^ 3 

nn/ iff ^(wim) n^.^(Kjf^-i) 



2iV 

X 



n n n (-^^ + ^j-* + ''^) n (-^i-* + + 

s<t t<s 



s<t s<t s>t 

n n ^("i.'' ~ "i:*)) n shQ(ao,s-ai,t)j J| sli Q(ai,t - ao,s) j . 

jr=0,l a,teAj ^ ^ sG^0''£-4i ^ ^ sg^g.te^l ^ ^ 



(5.6) 



Here we have set the kernel function $({a}|{/3}) by 



2N 

= nn H "Pif^^ - ^J,t) n ^(ao,s-ai,t)0(ao,s 
s=i j=o,it£Aj seAo,t€Ai 

^ n n - aj,t)(t>{aj,s - aj^t), (5.7) 



j=0,l s,teAj 

s<t 



where 



\37r 3y V Svr 3^ 
In the limiting case, our integral representation reproduce those of the correlation function at 
critical point q = —1 Q. 

Recently Kitanine, Maillet, Terras derived integral representation of the correlation 
function for the six-vertex model {A^^^ symmetry) by means of the quantum inverse scatter- 
ing method. Their method is avairable for both massive and massless regime. To construct 
A^^ analogue of the paper |^] and reproduce the integral representation ( ^.37 ) is interesting 
problem. 
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